Measurement of 03 with Dalitz Plot Analysis of B^ — > D^*'>K^ Decay 

A. Poluektov,! K. Abe/ T. Abe/ H. Aihara/^ y. Asano/^ T. Aushev," T. Aziz/^ S. Bahinipati/ A. M. Bakich/^ 
I. Bedny/ U. Bitenc/^ l. Bizjak/^ S. Blyth/^ A. Bondar/ A. Bozek/^ M. Bracko/^'^^ j^ Brodzicka/^ 
T. E. Browder,6 P. Chang/" Y. Chao/" B. G. Cheon/ R. Chistov," S.-K. Choi/ Y. Choi/^ A. Chuvikov/^ 
M. Danilov/i M. Dash/^ L. Y. Dong/ S. Eidelman/ V. Eiges/^ Y. Enari/" S. Fratina/^ N. Gabyshev/ 
A. Garmash/2 T. Gershon/ G. Gokhroo/^ B. Golob/^' ^^ j^ Guo/^ J. Haba/ T. Hara/^ K. Hayasaka/" 
H. Hayashii/i M. Hazumi/ T. Higuchi/ L. Hinz/^ T. Hokuue/° Y. Hoshi/o W.-S. Hou/'^ Y. B. Hsiung/"' * 
A. Imoto/i K. Inami/o A. Ishikawa/ R. Itch/ H. Iwasaki/ M. Iwasaki/^ j. H. Kang/^ j. S. Kang/" 
S. U. Kataoka/i N. Katayama/ H. Kawai/ T. Kawasaki/^ H. R. Khan/^ H. J. Kim/^ J. H. Kim/'^ S. K. Kim/^ 
K. Kmoshita/ P. Koppenburg/ S. Korpar/^'^^ P. Krizan/'''' ^^ P. Krokovny/ A. Kuzmm/ Y.-J. Kwon/® 
S. H. Lee/5 rp_ Lggi^k/s J. Li/" S.-W. Lin/" J. MacNaughton/° F. Mandl/° D. Marlow/^ T. Matsumoto/" 
. . . , A. Matyja/5 W. MitaroflF/° K. Miyabayashi/^ H. Miyata/"^ R. Mizuk/^ T. Mori/^ T. Nagamine/^ Y. Nagasaka/ 
■^ ■ T. Nakadaira/2 M. Nakao/ H. Nakazawa/ Z. Natkaniee/^ S. Nishida/ O. Nitoh/^ T. Nozaki/ S. Ogawa/^ 

§ ! T. Ohshima/'' T. Okabe/" S. Okuno/^ S. L. Olsen/ W. Ostrowicz/5 H. Ozaki/ P. Pakhlov/i N. Parslow/^ 

^ ; L. E. Piilonen/^ N. Root/ M. Rozanska/^ H. Sagawa/ S. Saitoh/ Y. Sakai/ T. R. Sarangi/ N. Sato/° 

,-H ■ O. Schneider/^ J. Schiimann/" C. Schwanda/" A. J. Schwartz/ S. Semenov/^ K. Senyo/° R. Seuster/ 

^ ; M. E. Sevior/^ H. Shibuya/^ B. Shwartz/ V. Sidorov/ A. Somov/ N. Soni/o R. Stamen/ S. Stanic/^' t 

^1 ■ M. Staric/2 A. Sugiyama/^ K. Sumisawa/^ T. Sumiyoshi/" S. Suzuki/^ O. Tajima/^ F. Takasaki/ K. Tamai/ 

! N. Tamura/"^ M. Tanaka/ G. N. Taylor/^ y. Teramoto/^ K. Trabelsi/ T. Tsuboyama/ T. Tsukamoto/ 

CN ; S. Uehara/ T. Uglov/i K. Ueno/" Y. Unno/ S. Uno/ G. Varncr/ K. E. VarveU/^ S. ViUa/^ C. H. Wang/^ 

^ ■ M.-Z. Wang/" M. Watanabe/^ B. D. Yabsley/"^ Y. Yamada/ A. Yamaguchi/i Y. Yamashita/^ M. Yamauchi/ 

^ ; J. Ying/i Y. Yusa/i S. L. Zang/ J. Zhang/ Z. P. Zhang/" V. Zhihch/ D. Zontar/^'^^ ^^^^ y). Zlircheri*^ 

O ; (The Belle Collaboration) 

jS^ ' ' Budker Institute of Nuclear Physics, Novosibirsk 

^1 I ^Chiba University, Chiba 

^L. . '^ Chonnam National University, Kwangju 

■"^.^ ' ^University of Cincinnati, Cincinnati, Ohio 45221 

^ ' ^ Gyeongsang National University, Chinju 

Qj , ''University of Hawaii, Honolulu, Hawaii 96822 

f») ■ ''High Energy Accelerator Research Organization (KEK), Tsukuba 

flj ' Hiroshima Institute of Technology, Hiroshima 

(~| , ^Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 

• • ■ ^''institute of High Energy Physics, Vienna 

^. ' ^''institute for Theoretical and Experimental Physics, Moscow 

k> , ^^J. Stefan Institute, Ljubljana 

5_] ■ ^^ Kanagawa University, Yokohama 

C^ ' ''^ Korea University, Seoul 

"^ Kyungpook National University, Taegu 
Swiss Federal Institute of Technology of Lausanne, EPFL, Lausanne 
'^'^ University of Ljubljana, Ljubljana 

''^University of Maribor, Maribor 
''"^ University of Melbourne, Victoria 

''^" Nagoya University, Nagoya 

^^ Nara Women's University, Nara 

^^ National Kaohsiung Normal University, Kaohsiung 

^'^ National United University, Miao Li 

''''^Department of Physics, National Taiwan University, Taipei 

^^H. Niewodniczanski Institute of Nuclear Physics, Krakow 

'''^Nihon Dental College, Niigata 

'^'Niigata University, Niigata 
^'^ Osaka City University, Osaka 

^'^ Osaka University, Osaka 

^'^Panjab University, Chandigarh 

''^ Peking University, Beijing 

'^'^ Princeton University, Princeton, New Jersey 08545 

' ' Saga University, Saga 

*'* University of Science and Technology of China, Hefei 



^^Seoul National University, Seoul 

Sungkyunkwan University, Suwon 

^'''University of Sydney, Sydney NSW 

^^Tata Institute of Fundamental Research, Bombay 

'^'^ Toho University, Funabashi 

^'^Tohoku Gakuin University, Tagajo 

'* Tohoku University, Sendai 

* Department of Physics, University of Tokyo, Tokyo 

^''^ Tokyo Institute of Technology, Tokyo 

^^ Tokyo Metropolitan University, Tokyo 

^^ Tokyo University of Agriculture and Technology, Tokyo 

^''University of Tsukuba, Tsukuba 

'^Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24-061 

'''^ Yonsei University, Seoul 

We present a measurement of the unitarity triangle angle 03 using a Dalitz plot analysis of the 
three-body decay of the neutral D meson from the B -^ D^*'K process. The method employs 
the interference between Z)° and 1)0 to extract both the weak and strong phases. We apply this 
method to a 140 fb~^ data sample collected by the Belle experiment. The analysis uses the modes 
B -^ DK and B — > D* K , D* — > Dtt", where the neutral D meson decays into Ks-iT'^n~ . 
We obtain 146 signal candidates for B -^ DK and 39 candidates for B — > D*K . From 
a combined maximum likelihood fit to the B -^ DK and B — > D*K modes, we obtain 
03 — 77° ^jgo(stat) ± 13°(syst) ± 11° (model). The corresponding two standard deviation interval 
is 26° < 03 < 126°. 

PACS numbers: 13.25.Hw, 14.40.Nd 

INTRODUCTION 

Determinations of the Cabbibo-Kobayashi-Maskawa (CKM) [1] matrix elements provide important checks on the 
consistency of the Standard Model and ways to search for new physics. The possibility of observing direct CP violation 
in _B -^ DK decays was first discussed by I. Bigi and A. Sanda [2]. Since then, various methods using CP violation in 
B -^ DK decays have been proposed [3-6] to measure the unitarity triangle angle 03 . These methods arc based on 
two key observations: neutral D° and -D" mesons can decay to a common final state, and the decay B^ -^ D^*^ K^ 
can produce neutral D mesons of both flavors via h -^ cus (Fig. la) and h —^ ucs (Fig. lb) transitions, with a relative 
phase 0+ between the two interfering amplitudes that is the sum, S + 4>3, of strong and weak interaction phases. For 
the charge conjugate mode, the relative phase is 0_ = (5 — 03, so both phases can be extracted from measurements of 
such charge conjugate B decay modes. However, the use of branching fractions alone requires additional information 
to obtain 03. This is provided either by determining the branching fractions of decays to flavour eigenstates (GLW 
method [3]) or by using different neutral D final states (ADS method [6]). 

A Dalitz plot analysis of a three-body final state of the D meson allows one to obtain all the information required 
for determination of 03 in a single decay mode. The use of a Dalitz plot analysis for the extraction of 03 was 
first discussed by D. Atwood, I. Dunietz and A. Soni in application to the ADS method [6]. This technique uses the 
interference of Cabibbo-favorcd Z?° -^ iir+7r~7r° and doubly Cabibbo-supprcssed D^ — > if +7r~7r° decays. However, the 
small rate for the doubly Cabibbo-suppressed decay limits the experimental applicability of this technique. Recently, 
three body final states common to 75° and D*^, such as Xstt+tt" [7], were suggested as being more promising, since 
both interfering amplitudes arc Cabibbo-favorcd in this case. This technique appears to have a higher statistical 
precision than methods based on branching fraction measurements. The statistical accuracy of the 03 extraction can 
be improved by adding the excited states of D and K to the analysis [8] . 

In the Wolfenstein parameterization of the CKM matrix elements, the amplitudes of the two diagrams that 
contribute to the decay i?+ -^ DK+ are given by Mi ~ Vcb^us ^ A\^ (for the D^K^ final state) and 
A/2 ^ ^ub^cs ^ A\'^{p + irj) (for D^K^). The annihilation diagram also contributes to A/2; but, since the weak 
coefficients are the same, this effectively leads to a redefinition of the strong phase. The two amplitudes A/i and 
A/2 interfere as the D^ and D^ mesons decay into the same final state /CgTr+Tr"; we denote the admixed state as D. 
Assuming no CP asymmetry in neutral D decays, the amplitude of the B^ decay can be written as 

M+ = /(m^,m?.) -hre*'^^+'*7("i-,"4), (1) 



i_(_ aiiu iit_ are the squared invariant masses of the Ks'k'^ and Kstt~ combinations, respectively, and 
f{m+,m-) is the complex amplitude for the decay -D° -^ iC^Tr+Tr". The absolute value of the ratio between the two 
interfering amplitudes, r, is given by the ratio l^^^bKsl/lK'^KtsI ^ 0.38 and the color suppression factor. The latter 
can be roughly estimated from the ratio of the color suppressed S" -^ D^K^ [9] and color allowed i3~ -^ D^K~ 



decays [10]: ^ B{B^ -^ D"KO)/B{B- -^ D"K-) = 0.35 ± 0.05. The amplitude ratio is therefore expected to be the 
product of these two factors, i. e. r ~ 0.13. However, other estimations of r exist, predicting the values as large as 
0.2 [11]. 

The corresponding amplitude for the charge conjugate B~ decay is 

Af_ = f{ml,ml) + re-''t''+'^ f{ml,ml_). (2) 

Once the functional form of / is fixed by a model for D'^ — > A'^tt+tt^, the D Dalitz distributions for 5+ and 
B^ decays can be fitted simultaneously using the above expressions for M+ and Af_, with r, 03, and 6 as free 
parameters. There are certain advantages of this technique: it is directly sensitive to the value of (f>3 and does not 
require additional assumptions on the values of r and S. Moreover, the value of r obtained in the fit can be useful for 
other (1)3 measurements. 

Reference [7] suggests a model-independent way for determining 03 via a binned Dalitz plot analysis. However, 
the application of this procedure would result in large statistical errors with our present data sample. Instead, we 
use a model-dependent approach based on an unbinned maximum likelihood fit to the D -^ Kstt'^tt" Dalitz plot 
distributions corresponding to 5+ and B^ data samples, thus making optimal use of our small number of signal 
events. The model of D^ -^ Kstt^tt^ decay in our approach is determined from a large sample of flavor-tagged 
D" — > KsTT^TT^ decays produced in continuum e^e^ annihilation. The drawback of this approach is that only the 
absolute value of the D° decay amplitude / can be determined directly; the complex form of / has to be based on 
model assumptions; these lead to potential model-dependent uncertainties in the determination of 03. Note, however, 
that the model uncertainties can be controlled in the future using data from r-charm factories. CP tagged neutral D 
mesons can be produced in the decay of the "0(3770) resonance, and these can be used to obtain information about the 
complex phase of the amplitude /, which is precisely the information required for a model-independent measurement 
of 03. 

The method used here has two possible two- fold ambiguities in the determination of the pair of parameters (03, S). 
The first one is a shift (03, (5) -^ (03 -|- tt, (5 -I- tt). The measured phases 9+ = S + (j)^ and 0_ = (5 — 03 do not change 
under this transformation. Another ambiguity is the inversion of sign (03,5) — > (— 03,— 15) with the simultaneous 
complex conjugation of the I?° decay amplitude /. This transformation does not change the observables, which are 
the squared absolute values of the amplitudes. However, if the D^ decay amplitude is approximated by a set of 
two-body amplitudes, the Breit-Wigner dependence fixes the sign of the imaginary part of the D° decay amplitude 
(the complex conjugate Breit-Wigner amplitude does not satisfy the causality requirement), and the second ambiguity 
is thus resolved. 

In a preliminary version of this analysis [12], we used only the B^ — > DK^ mode to constrain 03. The current 
measurement is based on two modes, B^ -^ DK^ and B^ —^ D*K^ {D* — > Dtt^). The statistical error evaluation 
is also improved compared to the previous measurement. 

EVENT SELECTION 

We use a 140 fb"-^ data sample collected by the Belle detector. The decays B^ -^ DK^ and B^ -^ D*K^, 
D* -^ Dtt^ are selected for the determination of 03; the decays B^ -^ Dtt^ , B^ -^ D*tt^ with D* -^ Dtt^ and 
B^{B^) -^ D*^7T^ with D*^ -^ Dn^ serve as control samples. We require the neutral D meson to decay to the 
Ksn^TT~ final state in all cases. We also select decays of D*^ -^ Dtt^ produced via the e+e~ — > cc continuum process 
as a high-statistics sample to determine the D" — + Kstt^tt^ decay amplitude. 

The Belle detector is described in detail elsewhere [13]. It is a large-solid-anglc magnetic spectrometer consisting 
of a three-layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC) for charged particle tracking 
and specific ionization measurement (dE/dx), an array of aerogel threshold Cerenkov counters (ACC), time-of- flight 
scintillation counters (TOF), and an array of 8736 CsI(Tl) crystals for electromagnetic calorimetry (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return located outside the 
coil is instrumented to detect Kl mesons and identify muons (KLM). 

Separation of kaons and pions is accomplished by combining the responses of the ACC and the TOF with the dE/dx 
measurement from the CDC to form a likelihood C{h) where ft, is a pion or a kaon. Charged particles are identified 



as pions or kaons using the likelihood ratio 7?.piD(ft-) = C{h)/{C{K) + C{-k)). 

Charged tracks are required to satisfy criteria based on the quality of the track fit and the distance from the 
interaction point in both longitudinal and transverse planes with respect to the beam axis. To reduce the low 
momentum combinatorial background we require each track to have a transverse momentum greater than 100 MeV/c. 
For charged kaon identification, we require the track to have TZpiYi{K) > 0.7. 

Photon candidates are required to have ECL energy greater than 30 MeV. Neutral pion candidates are formed from 
pairs of photons with invariant masses in the range 120 to 150 MeV/c^, or less than two standard deviations from the 
7r° mass. 

Neutral kaons are reconstructed from pairs of oppositely charged tracks without any pion PID requirement. We 
require the reconstructed vertex distance from the interaction point in the plane transverse to the beam axis to be 
more than 1 mm and the invariant mass A/^tt to satisfy IMtttt — Mks\ < 10 McV/c^, or less than four standard 
deviations from the nominal Ks mass. 

Selection of £>*=^ -^ Dtt"^ 

To determine the C decay model we use D* mesons produced via the e^e^ -^ cc continuum process. The flavor 
of the neutral D meson is tagged by the charge of the slow pion (which we denote as ng) in the decay D*^ -^ Dnf . 

To select neutral D candidates we require the invariant mass of the KsTr~^7T~ system to be within 9 MeV/c^ 
of the D° mass, Mdq. To select events originating from a D*^ decay we make a requirement on the difference 
AM ~ A'/^gjr+Tr-iTa ~ MxsTT+TT- of the invariant masses of the D*^ and the neutral D candidates: 144.6 McV/c^ < 
AM < 146.4 McV/c^. To suppress combinatorial background from BB events, we require the D* to have momentum 
in the centcr-of-mass (CM) frame greater than 2.7 GeV/c. 

The distributions of AAI and Af^^^+jr- for these events are shown in Fig. 2. The signal region bounds arc indicated 
with dashed lines. The resolutions of the selection variables are o'(AAf) = 0.38 MeV/c^ and fT(Afxs-n-+-n--) = 5.4 
MeV/c^. The number of events that pass all selection criteria is 104204. To obtain the number of background 
events in our sample we fit the AAf distribution. The background is parameterized with the function 6(AM) ^ 
(l/AA/)Y/(AA//m7r)^ — 1; the function describing the signal is a combination of two Gaussian peaks with the same 
mean value. The fit yields 100870 ± 840 signal events and 3210 ± 50 background events corresponding to a background 
fraction of 3.1%. 

Selection of B^ -^ DK^ 

The selection of B candidates is based on the CM energy difference AE = ^ E"; — i?bcam and the beam-constrained 
B meson mass A/bc = \/-^beam ~ (S-PO^j where i^bcam is the CM beam energy, and Ei and pi are the CM energies 
and momenta of the B candidate decay products. We select events with Mbc > 5.2 GeV/c^ and \AE\ < 0.2 GeV 
for the analysis. The requirements for signal candidates are 5.272 GeV/c^ < Mbc < 5.288 GeV/c^ and \AE\ < 0.022 
GeV. In addition, we make a requirement on the invariant mass of the neutral D candidate: IMkstttt — Muo\ < 11 
MeV/c2. 

To suppress background from e^e^ ^ 99 (<? = u,d,s,c) continuum events, we require | cos^thrl < 0.8, where ^thr 
is the angle between the thrust axis of the B candidate daughters and that of the rest of the event. For additional 
background rejection, we use a Fisher discriminant composed of 11 parameters [14]: the production angle of the B 
candidate, the angle of the B thrust axis relative to the beam axis and nine parameters representing the momentum 
fiow in the event relative to the B thrust axis in the CM frame. We apply a requirement on the Fisher discriminant 
that retains 90% of the signal and rejects 40% of the remaining continuum background. 

The AE and Afbc distributions for B^ -^ DK^ candidates are shown in Fig. 3. The peak in the AE distribution 
at AE = 50 MeV is due to B^ -^ Dtt^ decays, where the pion is misidentified as a kaon. The ratio of the number 
of events in the peak at AE = 50 MeV and in the signal peak is 0.54 ± 0.11, which is consistent with the ratio 
of B^ -^ DK^ and S± -^ Dtt^ branching fractions of 0.079 ± 0.009 ± 0.006 [10] and a 5% tt / K misidentification 
probability for our TIpyd{K) requirement. The B^ -^ DK^ selection efficiency (11%) is determined from a Monte 
Carlo (MC) simulation. The number of events passing all selection criteria is 146. The background fraction is 
determined from a binned fit to the AE distribution, in which the signal is represented by a Gaussian distribution 
with mean AE = 0, the B^ -^ Dn^ component is represented by a Gaussian distribution with mean AE = 50 MeV 
and the remaining background is modeled by a linear function. The contributions in the signal region are found to be 



112±12 signal events, 1.1±0.2 B^ -^ Dn^ events and 35±3 events in the linear background. The overall background 
fraction is 25 ± 4%. 

Selection of B± -^ D* K^ 

For the selection of B^ — > D*K^ events, in addition to the requirements described above, we require the mass dif- 
ference AM = M^s-K+Tv-Tvo -Mks^+^- of neutral D* and D candidates to satisfy 140 MeV/c^ < AM < 145 MeV/c^ 
Figure 4 shows the AE, Mbc and AM distributions for B^ -^ D*K^ candidates. The selection efficiency is 6.2%. 
The number of events satisfying the selection criteria is 39. The background fraction is determined in the same way 
as for B^ -^ DK^ events. The fit of the AE distribution yields 34 ±6 signal events, 4.4 ± 1.1 events corresponding to 
the linear background and 0.24 ± 0.08 B^ -^ D*Tr^ events in the signal region. The background fraction is 12 ± 4%. 

DETERMINATION OF D° -* Ksn+n- DECAY MODEL 

The amplitude / of the D^ -^ A'stt+tt^ decay is represented by a coherent sum of two-body decay amplitudes plus 
one non-resonant decay amplitude, 

Af 

f{ml,m'i) ^^aje'"'^ Aj{m\,m^_) +be'^ , (3) 

i=i 

where N is the total number of resonances, Aj{m?^,rn?_)^ qj and aj are the matrix element, amplitude and phase, 
respectively, of the j'-th resonance, and b and /3 are the amplitude and phase of the non-resonant component. The total 
phase and amplitude are arbitrary. To be consistent with a CLEO analysis [15], we have chosen the I?° — > Ksp mode 
to have unit amplitude and zero relative phase. The description of the matrix elements follows Ref. [16]. The matrix 
elements for the resonances are parameterized by Breit-Wigner shapes with D meson and intermediate resonance 
form factors and angular dependences taken into account. If we consider the decay of D^ into a particle C and a 
resonance r, with spin J, that subsequently decays into particles A and B, the expression for the matrix element is 

A = FoFr ""^ 



il/2 - Mlg - iMrTAB ' 



where Mr is the mass of the resonance, Mab is the invariant mass of the AB system, Fjj and Fr are the form factors 
of the D^ and the resonance, respectively, F^^ is the mass dependent width of the resonance, and sj accounts for the 
angular momentum of the resonance. The form factors are the Blatt-Weisskopf penetration factors [17]; both depend 
on the spin J of the intermediate resonance. We use the following expressions for the form factors: 



F= 1 



for J = 0, 



' 1 + R^P^r 



for J = 1, and 



/ 9 + 3i?2p2 + i?4p. 



9 + 3R^p\g + R'^p^ 



AB 



for J — 2. Here _R is a parameter that describes the radial size of the meson (either D*^ or resonance r), pab is the 
daughter particle momentum in the meson rest frame and Pr is its value when Mab = M^ (for Fjj the daughters are 
C and the resonance r, for Fr the daughters are A and B). The radial parameters we use are R ~ 5 GeV~^ for the 
D^ and i? = 1.5 GeV~^ for all intermediate resonances. The mass dependent width is given by 



2.7+1 



^--M't') m^'' 



where F^ is the width of the resonance. The angular term sj depends on the spin of the resonance. The expressions 
for scalar, vector and tensor states are: 

So = 1, 

Sl — ^^^AC ~ '^^BC '^ 



KB 



,,,2 ,,2 , (A/i,-M^)(Mi-Mi) 
S2 = 1 Mbc - ^"Ac H j;^ 



1 ( j,,r2 os/r2 oA.r2 , (^'^D ^-^c) \ /'n.r2 oA,r2 o a r 2 , (^^A 



- \^MXb ~ 2Af^ - 2M^ + - "^^^, "- j (^Mi^ - 2M1 - 2A/^ , ^^ 

For the I?° model we use a set of 15 two-body amplitudes. These include four Cabibbo-allowed amplitudes: 
ii'*(892)+7r-, A'f!;(1430)+7r-, i^|(1430)+7r- and A'*(1680)+7r-; doubly Cabibbo-suppressed partners for each of these 
states; and seven channels with Kg and a tttt resonance: Ksp, KsuJ, Ksfo{Q80), 1^5/2(1270), ivrs/o(1370), Kgcri 
and Ks<J2- The masses and Breit-Wigner widths of the scalars ci and 172 are left unconstrained, while the masses 
and widths of other resonances are taken to be the same as in the CLEO analysis [15]. In contrast to the CLEO 
analysis, we have introduced all doubly Cabibbo-suppressed amplitudes for flavor-specific decays (only 7^*(892)^7r+ 
was considered by CLEO) and two scalar states cti and (72 • The amplitude for cti describes the excess of events near 
the low TTTT invariant mass edge of the phase space. The resonance (T2 was introduced to describe a structure near 1.1 
GeV^/c'' in the m^^ distribution. This structure could be due to the decay /o(980) — > 7777 with rescattcring of 7777 to 
7r+7r~, which could distort the /o(980) — *■ tt+tt" amplitude for 77i7r7r near the 7777 production threshold. 

We use an unbinned maximum likelihood technique to fit the Dalitz plot distribution to the model described by 
Eq. 3. We minimize the inverse logarithm of the likelihood function in the form 



-21ogL = -2 



71 p 

\^ logj>(m^^^, 77T,_^^) — log / p{m^,'ni_)dm_^dm_ 



(4) 



D 

where i runs over all selected event candidates, and m'^ ^, m?_ ^ are measured Dalitz plot variables. The integral in 
the second term accounts for the overall normalization of the probability density. 
The Dalitz plot density is represented by 



p{m\,mi) ^ e{m\,mt) / |A/(m^ + ^^ 771^ + /i^)p cxp f ~ ^ ^^ ^ \ d^? + B{m\,m^_), 



(5) 



where M{m\, m?_) = f{m'^,mt) is the decay amplitude described by Eq. 3, e{m'^, m^) is the efficiency, B{m^, m?_) 
is the background density, amirn^^) is the resolution of the squared invariant mass m^^ of two pions (m^^ = 
M^ + M]^ + 2M^ — 777^ — 777^). The free parameters of the minimization are the amplitudes Qj and phases aj of the 
resonances (except for the Ksp component, for which the parameters are fixed), the amplitude b and phase /3 of the 
non-resonant component and the masses and widths of the cti and a2 scalars. 

The background density for D^ — > Kstt^tt" events is extracted from AM sidebands: AM < 142 MeV/c^ and 148 
MeV/c^ < AM < 150 MeV/c^. The background density is parameterized by a third-order polynomial in the variables 
777.^ and m^ to describe the purely combinatorial background, plus Dalitz plot densities for Z?" and D^ decays that 
correspond to events where a correctly reconstructed D^ is combined with a random slow pion. From the fit, we 
obtain the fractions of the background components: the purely combinatorial background is 43 ± 3%, combinations 
of Z?° with a pion of the correct charge account for 49 ± 3%, and the remaining 8 ± 1% is due to D'^'s combined with 
a pion of the wrong charge. The background fraction is fixed to 3.1% from the fit to the AM distribution. As a 
consistency check, we also perform a fit with the background fraction floated and obtain a value for the background 
fraction in agreement with the result from the AM fit. 

Our analysis is not sensitive to the absolute value of the reconstruction efficiency, but variations of the efficiency 
across the phase space can effect the fit result. The shape of the efficiency over the Dalitz plot is extracted from a 
MC simulation, where the D^ decays uniformly over the allowed phase space. The parameterization of the efficiency 



shape is a third-order polynomial in the variables m^ and m?_ , and symmetrical under interchange of tt"*" and 7r~ . 
The efficiency is nearly uniform over the central part of the Dalitz plot, but drops by between 5% and 13% at the 
edges of phase space. Finite momentum resolution has to be taken into account in the fit function since our model 
includes a narrow uj — > tt+tt" state. Although both m^_ and rn^ have finite resolution, we consider only the resolution 
of the 771^^ combination, since the other Dalitz plot projections do not contain any narrow structures. The resolution 
(Tm as a function of tti^tt is parameterized by a linear function and is extracted from MC simulation. The average 
ml^ resolution is 4.8 x lO'^ GeV^/c^. 

The fit results are given in Table I. The parameters of the a resonances obtained in the fit are: Mc^i = 539 ± 9 
MeV/c2, r,,^ = 453 ± 16 MeV/c^, M^,_ = 1048 ± 7 MeV/c^, and T^, = 109 ± 11 MeV/c^. The D" -^ Kstt+tt- 
Dalitz plot, as well as its projections with the fit results superimposed, are shown in Fig. 5. The large peak in the 
distribution corresponds to the dominant D^ -^ A'*(892)+7r~ mode. The minimum in the m^ distribution at 



m 



+ 

0.8 GeV^/c"* is due to destructive interference with the doubly Cabibbo suppressed I?^ -^ K*{892)~tt'^ amplitude. 
In the m^^ distribution, the I?° -^ Ksp contribution is visible around 0.5 GeV^/c"' with a steep edge on the upper 
side due to interference with D^ -^ Ks^- The minimum around 0.9 GeV^/c"' is due to the decay D° ^ A'5/o(980) 
interfering destructively with other modes. 

We obtain a larger amplitude for the non-resonant component compared to the CLEO analysis [15] (the fit fraction 
corresponding to the non- resonant component in our case is 24%) . The non-resonant component is found to be highly 
correlated with the amplitude for the ai resonance. A fit with the non-resonant amplitude fixed to zero yields a ai 
amplitude of 0.78 ± 0.05, while a fit without the ai yields a non-resonant amplitude of 4.66 ± 0.15. Therefore, we 
conclude that the large non-resonant fraction in our D^ model is due to a deficiency in our description of the ui state. 
We include this effect in the model uncertainty by performing additional fits to the B^ — » D^*^K^ data with D^ 
models with the non-resonant or cri amplitudes excluded. 

The unbinned likelihood technique does not provide a reliable criterion for the goodness of fit. To check the quality 
of the fit, wc make use of the large number of events in our sample and perform a binned x^ test by dividing the 
Dalitz plot into square regions 0.05 x 0.05 GcV^/c"'. The test yields x^ ~ 2121 for 1130 degrees of freedom. More 
detailed studies are required in order to understand the precise dynamics of D^ —^ Kstt'^tt^ decay. However, for the 
purpose of measuring (j)^ , we take the fit discrepancy into account in the model uncertainty. 

DALITZ PLOT ANALYSIS OF B± -* DK^ DECAY 

The Dalitz plots for D — > Kstt^tt^ , which contain information about CP violation in B decays, are shown in Figs. 6 
and 7 for B^ -^ DK^ and B^ — > D*K^, respectively. These distributions are fitted by minimizing the combined 
logarithmic likelihood function 

-21ogL = -21ogL_ -21ogL+, 

where L_(i+) are the likelihoods oi B^ (B^) data given by Eq. 4. The corresponding Dalitz plot densities p±(77z^, m?.) 
are given by Eq. 5 with decay amplitudes M± described by Eq. 1 (i?+ data) and Eq. 2 {B^ data). The D^ decay 
model / is fixed, and the free parameters of the fit are the amplitude ratio r and phases 03 and S. 

As in the study of the sample from continuum D* — > Dtt decays, the efficiency and the momentum resolution were 
extracted from the signal MC sample, where the neutral D meson decays according to phase space. The determination 
of the background contribution is described below. 

Backgrounds 

Five sources of background are considered in our analysis (see Table II). We determine the fraction and Dalitz 
plot shape for each component and use the results in the fit to the signal Dalitz plot. The largest contribution comes 
from two kinds of continuum events: random combination of tracks, and correctly reconstructed neutral D mesons 
combined with random kaons. These backgrounds are analyzed using an off-resonance data sample collected at an 
energy 60 MeV below the T(45') resonance in addition to a sample in which we make requirements on | cos^thrl and 
the Fisher discriminant that select continuum events and reject BB events. The continuum background fraction is 
22.1 ± 3.9% for B± -> DK"^ and 9.0 ± 3.6% for B^ ^ D*K^. The Dalitz plot shape of the continuum background is 
parameterized by a third-order polynomial in the variables m^ and Tm?_ for the combinatorial background component 
and a sum of D° and D*^ shapes for real neutral D mesons combined with random kaons. 



The background from BB events originates either from a B^ -^ D^*^K^ decay with some of the final state 
particles replaced by the decay products of the other B meson, or from other charged or neutral B decays (possibly 
with misidentified or lost particles). We subdivide the BB background into four categories: 

1. Decays other than B^ -^ D^*^ K^ and B^ -^ D^*\^, of which the dominant fraction comes from the decay of 
D^*"^ from one B meson, with some particles taken from the other B decay, constitute the largest part of the 
BB background. They are investigated with a generic MC sample. For this, we obtain background fractions of 
2.2 ± 0.2% for B± -^ DK^ and 2.1 ± 0.4% for B± -^ D*K^. The parameterization of this background includes 
a linear function in the variables rn?^ and ni^_ and a Gaussian peak in rn?_ . 

2. The process B -^ D^*>tt with a pion misidentified as a kaon is suppressed by the requirement on the K/tt 
identification probability and on the CM energy difference, AE. The fraction of this background is obtained by 
fitting the AE distribution; the corresponding Dalitz plot distribution is that of D" without the opposite flavor 
admixture. The fractions for this background are 1.0 ± 0.2% for B^ -^ DK^ and 0.6 ± 0.2% for B^ -> D* K^ . 

3. B^ -^ D^*'K events where one of the neutral D meson decay products is replaced by a random kaon or pion 
were studied using a MC data set where one of the charged B mesons from the T(4S') decays into the D^*'^K 
state. The corresponding background fraction is 0.4 ±0.1% for both B^ -^ DK^ and B^ — > D*K^ modes; the 
Dalitz plot shape is parameterized by a linear function in the variables m^ and ttj^ plus a D^ amplitude. 

4. Events in which a correctly reconstructed neutral D is combined with a random charged kaon are of importance, 
because half of the kaons have the wrong sign: such events will be misinterpreted as decays of D mesons of the 
opposite flavor, and thus introduce distortion in the most sensitive area of the Dalitz plot. In the MC sample, 
we flnd no events of this kind, which allows us to set an upper limit of 0.4% (at 95% CL) on the fraction for 
this contribution. 

Control sample fits 

To test the consistency of the fltting procedure, the same fitting procedure was applied to the B^ -^ D^*'>7r^ and 
B'^{B^) -^ D*^TT^ control samples as to the B^ -^ D'^*^K^ signal. For decays in which only one flavor D meson can 
contribute, the flt should return values of the amplitude ratio r consistent with zero. In the case of B^ — > D^*^!^^ 
a small amplitude ratio is expected [r ^ \VubV*^\ / \VcbV*^\ ~ 0.02). Deviations from these values can appear if the 
Dalitz plot distribution is not well described by the fit model. 

For the control sample fits, we consider i?+ and B~ data separately, to check for the absence of CP violation. The 
free parameters of the Dalitz plot fit are r± and 9±, where d± ~ 5 ± (ps (see Eqs. 1,2). 

The fit results for B^ ^ Dtt^ are r+ = 0.056 ± 0.028, 6+ = 237° ± 27° for B+ data and r_ = 0.068 ± 0.026, 
6- = 232° ± 22° for B^ data. It should be noted that since the value of r is positive definite, the error of this 
parameter does not serve as a good measure of the r = hypothesis. To demonstrate the deviation of the amplitude 
ratio r from zero, the real and imaginary parts of the complex amplitude ratio re'^ are more suitable. Figure 8 (a) 
shows the complex amplitude ratio constraints for the B~^ and B~ data separately. It can be seen from the plot that 
both amplitude ratios differ from the expected value by more than two standard deviations. This deviation is treated 
as a potential systematic effect. 

The other control samples, B^ -^ I)*7r± with D* decaying to I)7r°, and B°{B°) -^ D*±7r=F with 1}*^ -> Dt:^, do 
not show any significant deviation from r = 0. The results of the fit to the B^ -^ D*^^ sample (351 events) are 
r+ = 0.041 ±0.069, 6+ = 163° ±100°, r_ = 0.057 ±0.054, 9^ = 340° ±65° and are shown in Fig. 8 (b); the results of 
the fit to the B"(B") ^ D*^Tr^ sample (517 events) are r+ = 0.017 ± 0.070, 9+ = 278° ± 133°, r_ = 0.026 ± 0.050, 
9_ = 225° ± 99° and arc shown in Fig. 8 (c). 

RESULTS 

The results of the separate 5+ and B^ data fits are shown in Fig. 9. The plots show the constraints on the complex 
amplitude ratio re'^ for the B^ — » DK^ and B^ — » D*K^ samples. The fit technique is the same as the one used 
for the control samples. It can be seen that in both signal samples a significant non-zero value of r is observed. A 
difference between the phases 9+ and 0_ is also apparent in both the B^ -^ DK^ and B^ -^ D*K^ samples, which 
indicates a deviation of (j)^ from zero. 



A combined unbinned maximum likelihood fit to the B~^ and B samples with r, (f)^ and 5 as free parameters yields 
the following values: r = 0.31 ±0.11, (/)3 = 86° ±17°, S = 168° ± 17° for the B± -^ JDK^ sample and r = 0.34 ±0.14, 
(/)3 = 51° ± 25°, 6 = 302° ± 25° for the 5± -^ D*K^ sample. The errors quoted here are obtained from the likehhood 
fit. These errors are a good representation of the uncertainties for a Gaussian likelihood distribution, however in our 
case the distributions are highly non-Gaussian. In addition, the errors for the strong and weak phases depend on the 
values of the amplitude ratio r (e.g. for r = there is no sensitivity to the phases). A more reliable estimate of the 
statistical uncertainties is obtained using a large number of MC pseudo-experiments as discussed below. 

Estimation of model uncertainty 

The model used for the D*^ -^ A^^tt+tt" decay is one of the main sources of systematic error for our analysis. The 
model is a result of the fit to an experimental Dalitz plot, however, since the density of the plot is proportional to 
the absolute value squared of the decay amplitude, the phase 4>{'m'i^ , m?_) of the complex amplitude is not directly 
measured. The phase variations across the Dalitz plot are therefore the result of model assumptions and their 
uncertainties may affect the D Dalitz plot fit from B^ —> D^*'>K^. 

To estimate the effects of the model uncertainties, a MC simulation is used. Event samples are generated according 
to the Dalitz distribution described by the amplitude given by Eq. 1 with the resonance parameters extracted from our 
fit of continuum £)° data, but to fit this simulated plot different models for /(m4_,77i_) are used (see Table III). We 
scan the phases (jj^ and S in their physical regions and take the maximum deviations of the fit parameters ((Ar)„iax, 
(A(/f)3),nax, and (A(5)niax) as model uncertainty estimates. The values for (Ar)inax, (A(?!)3)niax and (A(S)max quoted in 
Table III are obtained with the value r = 0.13. For larger r values, the model uncertainty tends to be smaller, so our 
estimate of the model uncertainty is conservative. 

All the fit models are based on Breit-Wigner paramcterizations of resonances as in our default model. Since a 
Breit-Wigncr amplitude can only describe narrow resonances well, the usual technique to deal with broad states is to 
introduce Blatt-Weisskopf form factors for the D^ meson {Fd) and intermediate resonance (Fr) and a q^-dependence 
of the resonance width F. These quantities have substantial theoretical uncertainties and might introduce a large 
model error. We have therefore used a fit model without Blatt-Weisskopf form factors and with a constant resonance 
width to estimate such an error. We have also used a model containing only narrow resonances (A* (892), p, doubly 
Cabibbo-suppressed A*(892) and /°(980)) with the wide ones approximated by the flat non-resonant term. The 
study of the model errors is summarized in Table III. Our estimate of the systematic uncertainty on 03 is 11°. 

Estimation of systematic errors 

In addition to the model uncertainty, there are other potential sources of systematic error, such as uncertainties in 
the background Dalitz plot density, efficiency variations over the phase space and possible fit biases. These are listed 
in Table IV for the B^ -^ DK^ and B^ — > D*K^ modes separately. The component related to the background 
shape parameterization was estimated by extracting the background shape from the Md sidebands and by using 
a flat background distribution. The maximum deviation of the fit parameters from the "standard" background 
parameterization was assigned as the corresponding systematic error. The effect of the uncertainty in the background 
fraction was studied by varying the background fraction by one standard deviation. 

A potentially dangerous background is caused by events with a random kaon, half of which would not have the 
correct charge. We set a 0.4% upper limit on this kind of background at 95% confidence level based on MC simulation. 
The effect of this background on the fit results was studied using a MC procedure similar to that used for investigating 
the model uncertainty. The bias in the fit parameters corresponding to a 0.4% fraction is negligible (0.7° for (/>3) in 
comparison to the systematic error due to background shape. 

As mentioned above, the efhciency shape and momentum resolution were extracted from MC simulation. To 
estimate their contributions to the systematic error, we repeat the fit using a flat efficiency and a fit model that does 
not take the resolution into account, respectively. The biases due to the efficiency shape differ for B^ -^ DK^ and 
B^ -^ D*K^ samples, but since we expect the values of the efficiency systematics to be close for the two modes, we 
assign the maximum value of the bias as the corresponding systematic error. 

The non-zero amplitude ratio observed in the B — > Dtt control sample can be cither due to a statistical fluctuation 
or may indicate some systematic effect such as background structure or a deficiency of the D^ decay model. Since 
the source of this bias is unknown, wc conservatively treat it as an additional systematic effect. The corresponding 
bias of parameters is estimated in the following way. Suppose the parameters in the re*^ plane are biased by a value 
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Ar in a certain direction, then the maximum bias of the total phases would be equal to {/S.9)max ~ Ar/r ^ 11°. 
Since 5 = (6*+ + 0-)/2 and (f)^ = {9+ — 6^)/2, the maximum biases of the strong and weak phases would also be 
equal to 11°. The maximum bias of r for the simultaneous fit to both flavors in the case of 03 ^ tt would be 
Ar = ^r2 + (Ar)2 - r - 0.006. 



Evaluation of statistical error 



We use a frequentist technique to evaluate the statistical significance of the measurements. This method requires 
knowledge of the probability density function (PDF) of the fitted parameters as a function of the true parameters. 

To obtain this PDF, we employ a "toy" MC technique that uses a simplified MC simulation of the experiment which 
incorporates the same efficiencies, resolution and backgrounds as used in the fit to the experimental data. This MC is 
used to generate several hundred experiments for a given set of r, 9^ and 9_ values. For each simulated experiment, 
Dalitz plot distributions are generated with numbers of events that nearly equal to the numbers of events observed in 
the data — 70 events for each B flavor for B^ -^ DK^ and 20 events for each B flavor for B^ -^ D*K^ — and the 
simulated Dalitz distributions are subjected to the same fitting procedure that is applied to the data. This is repeated 
for different values of r, producing distributions of the fitted parameters that are used to produce a functional form 
of the PDFs for reconstructed values for any set of input parameters. 

We parameterize the PDF of a set of fitted parameters (r, (f>3,S), using the following model. We assume that as a 
result of the fit of a single Dalitz plot (either 5+ or B^ data), the errors of parameters i?e(r±e'^±) and Im{r±e^^^) 
are uncorrelated and have Gaussian distributions with equal RMS which we denote as a. The PDF of the parameters 
(r±, 9±) for the true parameters (f±, 9±) is thus written as 



d^P{r±,9±\f±,9±) 



1 



27rcr2 



exp 



{r± cos 6± — r cos 9±)^ + (r± sin 9± ~ r sin 9± 



i21 



2a2 



r±dr±d9±. 



To obtain the PDF for the parameters (r, 03 , S) we fix r 
6 - 4>3: 



and substitute the total phases with S - 



drd(f>3dS 



{r,(j)3,6\r,(j)3,S) 



d^P 
dr+d9^ 



■{r,S + (j)3\r,6 + (j)3) 



d^P 
'db~dBZ 



(7-,(5-03|r,5-03). 



and 



(6) 



In this model, there is only one free parameter a which is obtained from the unbinned maximum likelihood fit of the 
MC distribution of fitted parameters to Eq. 6. The values of a obtained for MC samples with different values of r are 
all equal within 4%, therefore we use a constant value a = 0.14 for the PDF of B^ -^ DK^ data. For the fit of the 
sample of neutral D mesons from B^ -^ D*K^ the value of a is equal to 0.22. 

After the PDF of the fitted parameters is obtained, the confidence level a for each set of true parameters (f , 03 , 5) 
is defined as 



a(r,03,(5) 



d^P 

Q drd(j)3dS 



[r, (p3, 5\r, (j>3, S)drd(l)3d5, 



where the corresponding confidence region Q is given by the condition 



d^P 

drd(f>3dS 



{r,(j>3,S\r,(j>3,S) > 



d^P 

drd(j)3dS 



(O.31,86°,168°|f,03,5), 



i.e. it includes all points in the fit parameter space for which the PDF is larger than that at the point 0.31, 86°, 168°, 
corresponding to the fit result. For B^ -^ D*K^ these values are replaced by 0.34, 51° and 302° for r, 03 and S, 
respectively. 

The confidence regions for the pairs of parameters (03,(5) and (03,?") are shown in Fig. 10 (B^ -^ DK^ mode) 
and Fig. 11 {B^ —> D*K^ mode). These plots are the projections of the corresponding confidence regions in the 
three-dimensional parameter space. We show the 20%, 74% and 97% confidence level regions, which correspond to 
one, two. and three standard deviations for a three-dimensional Gaussian distribution. The 20% confidence region, 
which corresponds to one standard deviation, yields the following results for the fit parameters: r ~ 0.26^Q'it4, 
03 = 86° ± 23°, (5 = 168° ± 23° for B± -^ DK^ data and r = 0.20t5^; J?, 03 = 51 ± 46°, (5 = 302 ± 46° for B± -> b*K^ 
data. The central values presented are obtained by maximizing the fit parameters' PDF. This technique accounts for 
the parameter biases introduced by the fit procedure. We find that 03 and 5 are unbiased in both the B^ -^ DK^ 
and B^ -^ D*K^ cases, while the central value of r is biased by the fit procedure due to its positive definitcness. 
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The values of the amplitude ratio r obtained are larger than our initial estimate (r = 0.13), though they agree 
within the statistical error. These values are also consistent with the recent measurement by BABAR collaboration 
[18], which set up an upper limit r < 0.22 at 90% CL for the B"^ -^ DK^ mode. 

In the frcqucntist approach, the significance of the CP violation is evaluated by finding the confidence level for the 
most probable CP conserving point, i.e. the point with r = or 03 = 0, for which the confidence level Oi{f^ 03, 5) is 
minimal. This procedure gives a = 97% for the point r = 0.03, 03 = 0, (5 = 168° for the B^ -^ DK^ sample. The 
same procedure applied to B^ -^ D*K^ sample gives a CP violation significance of 23% (for the point r = 0.10, 
03 = 0, (5 = 302°). 

Combined 03 measurement using B — > DK and B — > D* K samples 

The two events samples, B^ — > DK^ and B^ -^ D*K^^ are combined in order to obtain a more accurate 
measurement of 03. The technique we use to obtain the combined measurement is also based on a frequentist 
approach. Since in general the values of the amplitude ratio r and strong phase 6 can differ for the two modes, we 
have five true parameters (03, fi, f2, ^i and ^2, where the indices 1 and 2 correspond to B^ -^ DK^ and B^ -^ D*K^ 
modes, respectively) and six reconstructed parameters (r, 03 and 5 for each of the two modes). Since in this case 
the physical range of the parameters is smaller than the range of the reconstructed parameters (03 values have to be 
equal for the two modes), the Feldman-Cousins [19] approach is used. 

The PDF for the reconstructed parameters is written as 

T^(^'/^) "" , ^r'^X^ (^1' (<^3)l, (^ll?-'!, 03, ^"l) , ^7^Jr^ (^2, (03)2, '^2|r2, 03, ^"2), 

ax aracp^ad aracp^ad 

where x = {dri, d(03)i, d6i, dr2, d(03)2, dS2) is a vector of the reconstructed parameters, and fi = (03, fi, f2. Si, 62) is 
a vector of the true parameters. 

The confidence level a for a vector of true parameters /i is defined as 

Q^(/^) ^ / —;—{x\fi)dx. 
Jn dx 

The confidence region £7 is given by the Feldman-Cousins likelihood ratio ordering: 

dP . , IdP . . ,, dP , , IdP . , ,, 

dx I dx dx I dx 

Here purest ix) is defined as the vector of true parameters that maximizes the PDF for a given set x of reconstructed 
parameters. 

The vector of the central values of the true parameters is given by a = and equals fitestixo)- The corresponding 
central value of 03 is 77°. The one standard deviation interval for 03 (which corresponds to the 3.7% confidence level 
for the case of a five-dimensional Gaussian distribution) is 03 = 77° ^]^go ; the two standard deviation (or 45% CL for 
a five-dimensional distribution) interval is 39° < 03 < 112°. These intervals include only the statistical error. The 
statistical significance of CP violation for the combined measurement is 95%. 

Since the B — > DK contribution dominates in the combined measurement, we use its value of the systematic 
uncertainty, which is 13°, as an estimate of the systematic uncertainty in the combined 03 measurement. The model 
uncertainty for the two modes is the same and amounts to 11°. The two standard deviation interval including the 
systematic and model uncertainties is 26° < 03 < 126°. 

CONCLUSION 

We report results of a measurement of the unitarity triangle angle 03 that uses a new method based on a Dalitz plot 
analysis of the three-body D^ decay in the process B^ -^ D^*^K^. The first measurement of 03 using this technique 
was performed based on 140 fb~^ data sample collected by the Belle detector. From the combination of B^ -^ DK^ 
and B^ -^ D* K^ modes, wc obtain the value of 03 = 77° ^j^go ± 13° ± 11°. The first error is statistical, the second 
is experimental systematics and the third is model uncertainty. The two standard deviation interval (including model 
and systematic uncertainties) is 26° < 03 < 126°. The statistical significance of CP violation for the combined 
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measurement is 95%. The method ahows us to obtain a value of the D^-D^ amphtude ratio r, which can be used in 
other 03 measurements. We obtain r = 0.26l!5;i^±0.03±0.04 for the B^ -^ DK^ mode and r = 0.20lg;i[?±0.02±0.04 
for the S± -^ D*K^ mode. 

The method has a number of advantages over the other ways to measure 03 [3]-[6]. It is directly sensitive to the 
value of 03 and has only the two- fold discrete ambiguity (03 + tt, (5 + tt). It does not involve branching fraction 
measurements and, therefore, the influence of the detector systematics is minimal. The statistical sensitivity of this 
technique is also superior in the presence of background since an interference term is measured. 
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Kstt'^'k decay. Errors are statistical only. The results for the ai, (T2 masses and widths are 



Intermediate state 


Amplitude 


Phase (°) 


K*{892)+TT- 


1.656 ±0.012 


137.6 ±0.6 


K*{892)-TT+ 


(14.9 ±0.7) X 10"^ 


325.2 ± 2.2 


A'o*(1430) + 7r- 


1.96 ±0.04 


357.3 ±1.5 


A'o*(1430)-7r+ 


0.30 ± 0.05 


128 ±8 


A'2*(1430) + 7r- 


1.32 ±0.03 


313.5 ±1.8 


7^2* (1430) -7r+ 


0.21 ± 0.03 


281 ±9 


ii'*(1680) + 7r- 


2.56 ± 0.22 


70 ±6 


A'*(1680)-7r+ 


1.02 ±0.2 


103 ±11 


Ksp° 


1.0 (fixed) 


(fixed) 


Ksuj 


(33.0 ±1.3) X 10"^ 


114.3 ±2.3 


KsM980) 


0.405 ± 0.008 


212.9 ±2.3 


KsMmo) 


0.82 ±0.10 


308 ±8 


Ks f2il270) 


1.35 ±0.06 


352 ±3 


Ksai 


1.66 ±0.11 


218 ±4 


Ks (T2 


0.31 ± 0.05 


236 ±11 


non-resonant 


6.1 ±0.3 


146 ±3 



TABLE XL Fractions of different background sources. 



Background source 


B± -^ DK^ 


B± -» D*K^ 


qq combinatorial 


22.1 ±3.9% 


9.0 ± 3.6% 


BB events other than B^ ~> D^'''>K^/%^ 


2.2 ± 0.2% 


2.1 ±0.4% 


B± _> £)(*)7r±with K/-R misID 


1.0 ± 0.2% 


0.6 ± 0.2% 


Combinatorics in D^ decay 


0.4 ±0.1% 


0.4 ±0.1% 


Combinatorial kaon in B* — > D'^*^ K^ decay 


<0.4% (95% CL) 


<0.4% (95% CL) 


Total 


25 ± 4% 


12 ± 4% 



TABLE IV: Contributions to the experimental systematic error. 



Source 


B± ^ DK^ 


B± -^b-K^ 


Ar 


A<^3 n 


A5(°) 


Ar 


A<^3 n 


Asn 


Background shape 
Background fraction 
Efficiency shape 
Momentum resolution 
Control sample bias 


0.017 
0.025 
0.004 
0.010 
0.006 


4.7 
0.1 
3.5 
2.5 
11 


2.3 
0.6 
1.2 
0.6 
11 


0.016 
0.015 
0.002 
0.010 
0.006 


1.5 
0.6 
3.5 
2.5 
11 


2.6 
0.9 
1.2 
0.6 
11 


Total 


0.032 


13 


11 


0.024 


12 


11 



TABLE III: Estimation of model uncertainty. 



Fit model 


(Ar)max 


(A,^3)max (°) 


(A<5)^ax (°) 


Fr=FD = l 

r(g^) = Const 
Narrow resonances plus non-resonant term 


0.01 
0.02 
0.03 


3.1 

4.7 
9.9 


3.3 
9.0 
18.2 


Total 


0.04 


11 
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FIG. 1: Feynman diagrams of (a) dominant _B+ — > D'^K'^ and (b) suppressed _B+ -^ D'^K'^ decays 
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FIG. 2: (a) AM and (b) Md distributions for the D* -^ DyrJ^candidates. Dashed lines show the signal region. The histogram 
shows the data; the smooth curve in (a) is the fit result. 
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FIG. 3: (a) AE and (b) Afbc distributions for the B^ — > DiiT^ candidates. Dashed lines show the signal region. The histogram 
shows the data; the smooth curve in (a) is the fit result. 
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FIG. 4: (a) AE, (b) A4"bc and (c) AM distributions for the B -^ D*ii' candidates. Dashed lines show the signal region. The 
histogram shows the data; the smooth curve in (a) is the fit result. 
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FIG. 5: (a) 77i+, (b) m_, (c) jn^^ distributions and (d) Dalitz plot for the D^^ -^ A'stt^tt decay from the D* -^ D-Kg process. 
The points with error bars show the data, the smooth curve is the fit result. 
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FIG. 6: Dalitz plots of £» ^ T^Tstt+tt- decay from (a) B+ -> DK+ and (b) B" -^ DR- 
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FIG. 7: Dalitz plots of £) -» T^sTr+Tr" decay from (a) B+ -^ D*K+ and (b) B' -» D-R- 
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FIG. 8: Constraint plots of the complex amplitude ratio re^* for (a) B± -» D-k^ , (b) B± -^ D'-k^ and (c) BO{B°) -^ D'^-k"^ 
decays. Contours indicate integer multiples of the standard deviation. Dotted contours are from B^{B°) data, dashed contours 
are from B^IB*^) data. 
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FIG. 9: Constraint plots of the complex amplitude ratio re' for (a) B — > DK and (b) B — > D* K decays. Contours 
indicate integer multiples of the standard deviation. Dotted contours are from B'^ data, dashed contours are from B^ data. 
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FIG. 10; Confidence regions for the pairs of parameters (a) (r, (/is) and (b) {4>3,5) for tiie B — > DK sample. 
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FIG. 11: Confidence regions for the pairs of parameters (a) (r, ^s) and (b) {(j>s,5) for the B — > D* K sample. 



